Double exchange-driven spin pairing at the (001) surface of manganites 
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The (001) surface of Lai-^Ca^MnOa system in various magnetic orderings is studied by first 
principle calculations. A general occurrence is that z 2 dangling bond charge - which is "invisible" 
in the formal valence picture - is promoted to the bulk gap/Fermi level region. This drives a 
double-exchange-like process that serves to align the surface Mn spin with its subsurface neighbor, 
regardless of the bulk magnetic order. For heavy doping, the locally "ferromagnetic" coupling is 
very strong and the moment enhanced by as much as 30% over the bulk value. 
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Although most efforts on the colossal magnetoresis- 
tance (CMR) materials typified by the Lai-^Ca^MnOs 
(LCMO) system are still concentrated on bulk properties, 
growing interest is being shown in the surface behavior 
. Knowledge of surface properties is essential not 
only to develop a perovskite manganite-based technol- 
ogy but also to determine fundamental phenomena and 
mechanisms of magnetoelectronic behavior. Indeed, the 
CMR effect occurs at high temperature, around the mag- 
netic ordering temperature, and a magnetic field of sev- 
eral Tesla is required to suppress the thermal magnetic 
disorder and produce the change in resistivity. Since high 
magnetic fields are generally unavailable in applications, 
alternative ways to trigger large low-field MR were con- 
sidered, such as with trilayer junctions Q and polycrys- 
talline samples ^]. The junctions are epitaxially grown 
along the [001] direction, and are made of a central in- 
sulating thin film of SrTiC>3 (the barrier), sandwiched 
by two metallic layers of Lao.67Sro.33Mn03 (LSMO). Ap- 
plying a low magnetic field, the tunneling conductivity 
can be switched by inducing a parallel (switch on) or 
anti-parallel (switch off) spin-orientation in the two elec- 
trodes. Taking advantage of their half metallicity gives a 
very large tunneling MR (TMR). 

Large low-field intergrain MR (IMR) M over a large 
temperature range has been observed in polycrystalline 
samples of LSMO, Cr0 2 , § and the double per- 
ovskite systems Sr2Fe(Mo,Re)Og, P,|lO| all of which are 
expected to be half metallic magnets. Magnetotunncling 
across grain boundaries, in which the relative orientation 
of the magnetization of neighboring grains is manuipu- 
latcd by an applied field, is believed to be the mechanism. 
In the IMR process, which may be the most promising 
for MR applications, there is mounting evidence that the 
state of the surfaces of the grains is important in the 
intergrain tunneling process. JHI For TMR it has long 
been clear that tunneling characteristics are strongly in- 
fluenced, perhaps even dominated, by the electronic and 



magnetic structure at the interface, and for IMR surface 
states have been suggested to play the central role. 

In the few experimental works present in the literature 
intrinsic difficulties have been reported in the process of 
obtaining clean, bulk-truncated surfaces, due to surface 
segregation that occurs during growth at high temper- 
ature j|, and strain effects induced by film-substrate 
mismatch g. Structural and electronic properties of 
the low-index surfaces (including the possibility of re- 
constructions) are still unknown, in spite of their impor- 
tance in establishing the half metallic nature of the CMR 
materials using photoelectron emission. (I]] However, ad- 
vancements in epitaxial growth and surface uniformity 
are being reported, jll| so a first fundamental step to- 
wards describing real surfaces consists in understanding 
how the intrinsic properties of the ideal unreconstructed 
surfaces differ from the respective bulk properties, i.e. 
how the bulk truncation in itself modifies the physics of 
this compounds. First- principle calculations are ideally 
suited for this aim, and in this paper we focus on sur- 
face spin ordering of the Mn-terminated (001) surface 
of Lai-zCa-rMnOa. Although this system shows an ex- 
tremely rich variety of magnetic phases for different level 
of doping jl2| , we identify a general, robust mechanism 
that should dominate the surface spin order for any dop- 
ing level x. In the regime of heavy doping x ~ 1/2, the 
surface-to-subsurface magnetic coupling is much stronger 
than in bulk, and the surface moment is enhanced by as 
much as 30%. 

Based on the growing understanding of the double ex- 
change (DEX) process in bulk manganites, |l^] it can be 
expected that the surface spin alignment will be strongly 
dependent on the Mn e g occupation. At the (001) sur- 
face, however, the e g degeneracy is broken: the x 2 — y 2 
orbital remains very strongly dpa hybridized with neigh- 
boring (in surface layer) O ions, but the z 2 orbital is left 
"dangling." The implications of this symmetry breaking 
were first glimpsed in the simple, undoped x=\ member 
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CaMn03, which has G-type AFM bulk ordering due to 
standard AFM superexchange between filled t2 g shells. 
These t2 g shells contain the nominal three electrons as- 
signed to Mn 4+ in the formal valence picture, with the 
e g formally empty. As it has been pointed out in other 
contexts, Q, the amount of actual d charge is not at all 
identical to the formal d n charge. For bulk phenomena 
however, this idealization usually gives a reliable broad 
picture of general behavior, including spin, charge, and 
orbital order. For the CaMn03 (001) surface, however, 
G-type spin order does not survive at the surface. In- 
stead, a flip of all spins in the surface layer occurs, driven 
by the appearance of Mn z 2 charge 0] that drives a dou- 
ble exchange (DEX) process that strives to align spins. 
This z 1 charge is present in the bulk, resulting from dpa 
mixing that draws a truly significant amount of e g charge 
into the "O 2p" bands: the 18 O p bands actually con- 
tain on the order of 1.5-2 electrons of Mn e g character. 
Some of this becomes a dangling bond band at the sur- 
face, lying in the bulk gap and driving the surface spin to 
flip. What we show in this paper is that this mechanism 
survives, and in fact is enhanced, as doping occurs. 

First-principles calculations have been performed 
within local-density approximation (LDA), employing a 
plane wave basis and Vanderbilt pseudopotentials p5| . 
A 30 Ryd cut-off energy and the exchange-correlation 
potential of Perdew and Zunger [|l6| was used. For the 
Lai-zCa^MnOa (001) surfaces we have studied, we used 
slab of nine atomic layers. For £=1/2, the stacking along 
z is made by alternating layers of La and Ca (see Fig. 
|l|) , retaining a mirror symmetry with respect to the cen- 
tral Mn layer. The artificial ordering of La and Ca lay- 
ers (which must be done somehow in a finite supercell) 
should not affect our conclusions, since they simply be- 
come ionized by contributing their valence electrons to 
the O and Mn bands. For the planar lattice constant of 
bulk Lai/2Cai/2Mn03 we obtain by energy minimization 
ao= 7.21 a.u, which is a reasonable value between the ex- 
perimental 7.35 a.u. for La 2 /3Ca 1 / 3 Mn03, and 7.05 a.u. 
for CaMn0 3 . We find the AFM phase favored by 15 
meV/Mn over the FM, which has a nearly half-metallic 
density of states. 

In Table | the calculated relative surface energies for 
x= 1/2 are reported. As can be seen from Fig. 1, there are 
two kinds of Mn-terminated surfaces, i.e. one with La in 
the subsurface layer (indicated as Mn-La), and one with 
Ca instead (Mn-Ca). Since (Table I) they give equivalent 
results, we will quote specific results for just one of them 
(La-Mn). Treating an inner region with the true spin, 
charge, and orbital order for x=l/2 (l7) is well beyond 
computational capabilities, but the behavior we identify 
is so robust that we expect it to be independent of bulk 
order. For each of the geometries of our nine layer sur- 
face, four spin arrangements on Mn are possible, labeled 
in Table Q by arrows on central (C), subsurface (SS) and 
surface (S) Mn, in this order. These are: surface and 



subsurface parallel, aligned or antialigned to the central 
layer, and surface and subsurface antiparallel, with the 
subsurface aligned or antialigned with the central layer. 
We find that S-SS spin alignment is strongly favored, 
with the most stable configuration having both the sur- 
face and subsurface layer spins antiparallel to the central 
Mn spin: (from third to top layer) 

The energies can be mapped onto an interlayer Ising 
model with three independent effective exchange con- 
stants: Js-sSi Jss-c and Js-c> the latter being a sec- 
ond neighbor coupling. Jss-c = -18 meV (AFM) is close 
to the exchange parameter obtained directly from the 
bulk calculation (Jbuik = -15 meV). The interaction be- 
tween Mn on first and third layers, Jc-S = 8 meV, is 
FM in sign and is related to the d z 2 the surface state 
discussed below. 

The most striking result for £=1/2 (Table |) is the 
positive, unusually large value of Js-SS = 53 meV, more 
than three times larger than, and opposite in sign to, the 
bulk AFM coupling. For comparison, for CaMn03 (x — 
1) the interlayer exchange constant at the surface was 29 
meV. (The bulk coupling for x=l is also different from 
:r=l/2, with Jbuik — -26 meV.) This large FM coupling, 
for both x=l and x=l/2, is the consequence of a very 
general characteristic of (001) surface formation. 

In Fig.|| the orbital-resolved density of states (DOS) of 
the Mn ions for the (001) surface in the most stable spin 
configuration (i.e. |J,|) is shown. Two surface Mn d z 2 
DOS peaks straddle the Fermi energy (Ep = 0), with a 
tail of occupied states that extends down to ~ -1.5 eV. 
These states are also visible on subsurface Mn and, for 
the occupied peak, on central Mn as well. Thus the sur- 
face has a deep surface d 2 2 resonance extending to the 
fifth layer below the surface. In the majority channel of 
the central ('bulk') Mn ion, d z 2 and d x 2_ y 2 orbitals con- 
tribute to the DOS at Ef, whereas in the minority chan- 
nel the only contribution comes from t2 g states. The d z 2 
dangling bond discussed above leads to the formation of 
the surface resonance as it does in CaMnOs. It is also ap- 
parent that the d^y bands are shifted upward in energy, 
so that the minority channel is depleted (i.e. the Mn at 
surface is fully polarized) and the dj.^ surface bands con- 
tribute to the DOS at Ef- The magnetic moment on the 
surface Mn (3.23 (1b) is 10% larger than on subsurface 
Mn (2.97 hb) and 30% larger than in 'bulk' central Mn 
(2.50 hb), but the total charge on Mn (~ 5.3 electrons 
using our definition) is nearly the same at surface and in 
the bulk. The increase of magnetization is due mostly 
to the d z 2 polarization with some contribution from the 
depletion of d xy states around E^- Also, a small intra- 
atomic charge readjustment occurs from d x 2_ y 2 and d xy 
to the polarized d z 2 orbital on the surface Mn ion. 

The resulting surface polarization can be visualized 
from the isosurfaces of the magnetization displayed in 
Fig||. Contributions coming from states that lie in the 
region within 0.3 eV below Ep are shown, i.e. the "core" 
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tig moments are not included in the subsurface and cen- 
tral Mn, whereas the surface Mn shows a combination of 
d z 2 and d xy spins; on subsurface Mn the d z 2 magnetiza- 
tion is mixed with some d x 2_ y 2. The double exchange 
effect between d z 2 orbitals on surface and subsurface Mn 
comes into play and leads to the strong FM coupling 
Js~ss=53 meV responsible for the spin alignment. On 
the central Mn with its antialigned spin, the magneti- 
zation is d 2 ,2_j,2-likc. (Unfortunately, present computa- 
tional limitations do not allow us to study an eleven layer 
slab, for which the central layer should be more bulklike.) 
Also evident in Fig. || is that a remarkably large fraction 
of this surface-induced magnetization lies in the O p w or- 
bitals of the surface layer. Polarization of the O ion in 
FM bulk environments in manganites has been empha- 
sized elsewhere. Q 

The change of the Mn d z 2 orbital from broad, strongly 
dpa hybridized in the bulk to an atomic-like, narrow in 
energy, surface state is a very specific feature of this (001) 
surface formation, and this surface dehybridization gen- 
erally should be described well by LDA. We suggest that 
this effect is strong enough to turn the AFM spin cou- 
pling into FM for any doping level. At least two argu- 
ments support this hypothesis. First, the spin-pairing 
occurs for the (001) surface of CaMn03 H| that should 
be the most unfavorable case, since in the bulk (nomi- 
nally) only the majority t2 9 orbitals are occupied, thus 
their AFM character is dominant. Nevertheless, the par- 
tially occupied d z 2 surface state reverses the magnetic 
coupling. Second, the very large change of exchange in- 
teraction parameter (from -15 meV in bulk to +53 meV 
at the surface) would overcome AFM bulk coupling even 
stronger than the one considered here. 

A crucial case is the x=0 member LaMnC>3, which is 
A-type AFM in the bulk. The spin-pairing argument 
applied to the surface parallel to the FM (001) layers 
predicts a spin-flip of the surface Mn layer. The AFM 
spin coupling along the z axis is robust and explained 
by a well estabilished picture: the in-plane FM coupling- 
is stabilized by the ordering of Mn e g orbitals, so that 
occupied d x 2 (d y 2) orbitals alternates with empty d x 2 
(d y 2) orbitals on neighboring Mn. Thus, all the e g -type 
charge fills in-plane orbitals, and the d z 2 orbitals are 
empty and higher in energy. As a consequence, the AFM 
interactions between neighboring t2 9 's dominates in the 
orthogonal direction. A realistic first-principle calcula- 
tion of the LaMnC>3 surface is beyond the possibility of 
detailed calculations, since it would require a 
lateral enlargement of the cell as well as additional thick- 
ness to treat the tilting of the MnC>6 octahedra |2(J , and 
the Jahn- Teller distortion at the surface would have to 
be determined. However, the formation of the d z 2 sur- 
face state within bulk LaMnC>3 gap seems to be beyond 
doubt, based on the behavior of the d z 2 dangling bond for 
x=l and £=1/2. The question is whether this would be 
able to overcome the t2 ff AFM contribution. In Ref. jni] 



it is shown that the t2 9 contribution for bulk LaMnC>3 
increases linearly in magnitude with the distortion be- 
tween in-plane and inter-planar lattice constants at fixed 
volume, i.e. the AFM coupling between (001) planes in- 
creases linearly by shortening the interplanar distance, 
likely due to the electrostatic repulsion that further de- 
pletes the d z 2 orbitals. The tig contribution to Jbuik has 
been calculated in Ref. |l9| as a function of the lattice 
distortion. It is in the range of ~ 20-30 meV, i.e. not 
large enough to overcome the value of Js-sSj thus we 
definitely expect the occurrence of a spin-flip process at 
the (001) LaMn0 3 surface. 

In Fig.[| we show the e g orbitals on surface and subsur- 
face layers, and indicate the expected filling and orbital 
ordering after the formation of the surface state. The 
orbitals are ordered in "FM" fashion both in-plane and 
orthogonally to the surface, as a result of the surface for- 
mation that fills the Mn surface d z 2 orbital no longer 
degenerate with the d z 2 orbital of the underliyng subsur- 
face Mn. 

To summarize, we have found that terminating the 
(001) surface of Lai-^Ca^MnOs surfaces with the Mn 
ion exposed, results in partial filling of the d z 2 orbital 
that drives a double-exchange-like ordering of the sur- 
face and subsurface layers of Mn ions. We have shown 
this effect explicitly for x=l/2 and (previously) for un- 
doped CaMnC>3. A comparison between these two cases 
indicates that it is stronger in doped systems. This result 
has important implications (1) for surface studies, where 
this effect tends to insure that surfaces of the CMR mate- 
rials (x ~l/3) will remain ferromagnetically aligned and 
half metallic as well, as supported by photocmission stud- 
ies, and (2) for the intergrain magnetoresistance effect, 
where the magnetic structure of the grain surfaces can 
strongly affect the device characteristics. This behavior, 
which is strongly related to band filling but much less 
dependent on ion size effects, should also hold for the 
Lai-^Sr^MnOs and Lai-^Ba^MnOa systems. 
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TABLE I. Energies (in meV) for different spin con- 
figurations on Mn atoms. Each of them is labeled by 
the three arrows indicating the spin direction of cen- 
tral-subsurface-surface Mn. Mn-La is the Mn-terminatcd 
(001) surface with La on second layer, Mn-Ca is that one 
with Ca on second layer. All energies refers to that one of the 
most stable arrangement, i.e. ][[. 
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FIG. 1. Structure of the (001) surface of tetragonal 
La^Ca^MnOs. Arrows indicate spin orientations for the 
stable HI ordering (see text), which we find for this com- 
pound and for CaMnOs. Parallel alignment of the surface 
and subsurface layers is expected to be true generally. Note 
that the surface oxygen ions are also polarized. 
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FIG. 2. Orbital-resolved Mn d DOS for the (001) Mn-La 
surface in the spin configuration fj.-l (i.e. AFM bulk and 
flipped spin at surface) . The panels refer to the 3 unequivalent 
Mn layers in the slab. 




FIG. 3. Isosurfaces of the valence magnetization for the 
(001) Mn-La surface in the spin configuration The mag- 
netization shown is due to states below, and within, 0.3 eV 
of Ei? (see the corresponding DOS in Fig.^). Dark and light 
isosurfaces are of same magnitude but opposite sign, i.e. they 
represent up and down spin densities, respectively. 
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FIG. 4. Orbital ordering at LaMn0 3 (001) surface: filling 
of d z 2 orbital (indicated by shading) at the surface produces 
FM spin coupling perpendicular to the surface. Planar orbital 
ordering in the subsurface and other buried layers leads to FM 
layers alternating in spin direction, except at the surface. 
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